Optical link synthesis system

ABSTRACT

A method of optimising the componentry of an optical data transmission link, the link comprising a series of optical components interconnected by a fiber network, the method comprising the steps of: (a) calculating an equivalent fiber representation for the link; (b) utilising an optimisation algorithm to eliminate redundant optical components from the link; (c) utilising an optimisation algorithm to determine the input channel power of the link.

FIELD OF THE INVENTION

[0001] The present invention relates to the creation of optical networks and, in particular, discloses various algorithms to optimise the placement of optical elements within a network.

BACKGROUND OF THE INVENTION

[0002] Recently, the construction of optical networks has become more and more significant for forming the basis for data transmission. Optical fibers have been found to have higher and higher levels of useable capacity and hence they are beginning to dominate the area of long haul data transmission.

[0003] The design of an optical fiber network is often an extremely complex process. The layout often includes very complex components that must be placed at certain geographical locations so as to meet overall performance requirements. Often, it is desirable to minimise the number of components in a network so that over engineering is avoided.

[0004] As the creation of an optical network can be an extremely complex and expensive process, there has risen a general need for the accurate simulation of optical networks utilizing software based systems. Simulations allow a network design to be “trialed” in a virtual environment, which is much faster and less expensive that constructing an experimental prototype, laboratory or field trial system. However, while simulation is a very cost effective method to verify an existing design, by itself simulation does not filly address the overall problem of system synthesis.

[0005] Currently, the initial design (or synthesis) of an optical transmission system is a complex task that is best carried out via the intensive attention of an expert human designer. It is highly desirable to be able to capture equivalent expertise within a software synthesis system, to enable optical transmission systems to be preferably designed by people lacking such high levels of expertise.

[0006] The process of synthesis involves taking an initial set of design constraints and objectives, and generating a system design using equipment selected from a specified range of suitable, commercially-available components. The design constraints typically include such factors as existing installed fiber plant and geographical facility locations that must be used. The design objectives typically consist of overall performance and cost requirements.

[0007] Such a software synthesis system allows for complex planning to be undertaken, including the testing of multiple alternatives without having to place or alter complex and expensive structures in the field.

SUMMARY OF THE INVENTION

[0008] It is an object of the present invention to provide for an improved form of optical network synthesis and planning.

[0009] In accordance with a first aspect of the present invention there is provided a method of calculating an equivalent fiber representation of a multiplicity of actual fibers, the method comprising the steps of successively combining the parameters of an initial pair of adjacent fibers to compute the parameters of a third fiber that has substantially the same transmission properties as the two original fibers, and wherein one of the initial fibers may itself be an equivalent fiber representation. Alternatively, the method may comprise the step of combining the parameters of a multiplicity of fibers, comprising an arbitrary combination of actual fibers and equivalent fiber representations, to compute the parameters of a single replacement fiber that has the same transmission properties as the multiplicity of fibers.

[0010] In accordance with a second aspect of the present invention, there is provided a method for optimizing the topology of an optical data transmission link, the link comprising a series of optical components interconnected by fiber sections, the method comprising the steps of:

[0011] a) calculating an equivalent fiber representation of the link;

[0012] b) eliminating redundant optical components from the link;

[0013] Preferably the method also contains the step of:

[0014] c) generating solutions utilizing components selected from a set of available devices;

[0015] Preferably the method also contains the step of:

[0016] d) utilising an optimisation procedure to determine the input channel power of the link.

[0017] The step (a) preferably can comprise the steps of: (a1) progressively joining segments together between facilities and determining whether an equivalent fiber representation for the joined segment falls within predetermined constraints; (a2) progressively joining segments together across facility locations and determining whether an equivalent fiber representation for the joined segment falls within predetermined constraints.

[0018] The step (b) can comprise the steps of: (b1) calculating a measure of the signal quality required at the receiver to satisfy the link performance criteria; (b2) calculating the maximum acceptable transmission power for said link; (b3) determining suitable configuration values for predetermined components in the link. The configuration values can, for example, include the amplification factor of an optical amplifier, and the total dispersion of a dispersion compensating module.

[0019] The step (b1) can further comprise the step of calculating a numerical measure of signal quality corresponding to a specified link performance. The measure of link performance can, for example, be the probability of transmission errors occurring on the link. The numerical measure of signal quality can be, for example, the Q-factor. The step (b1) can also comprise the step of adjusting this signal quality measure to account for additional sources of degradation in transmission. The additional sources of degradation can be, for example, nonlinear transmission processes such as four-wave mixing. The step (b1) can further comprise the step of calculating a simplified measure of signal quality at the receiver that is more easily optimised in the subsequent steps. The simplified measure of signal quality can be, for example, the optical signal-to-noise ratio (OSNR).

[0020] The step (b2) can further comprise the step of calculating the maximum transmission power that will not result in unacceptable signal degradation due to nonlinear transmission effects and additive noise. The nonlinear transmission effects considered can include, for example, four-wave mixing and stimulated Raman scattering. The additive noise considered can be, for example, amplified spontaneous emission noise introduced in optical amplifiers. The step (b2) can further comprise the step of verifying that the output power is valid, within the operating ranges of components in the link. This can include, for example, verifying that the transmitters are capable of generating the required output power.

[0021] The step (b3) preferably includes the step of determining the amount of residual dispersion required to optimise the signal quality at the receiver, considering the impact of nonlinear propagation effects. The nonlinear propagation effects can be, for example, the self phase modulation of each transmitted channel. The step (b3) preferably further includes the step of applying dispersion compensation to achieve the required residual dispersion by selecting from available dispersion compensating modules. The available dispersion compensating modules can include, for example, lengths of dispersion compensating fiber. The step (b3) can also include the step of compensating for the insertion loss of the dispersion compensating modules by including further optical amplifying units in the link.

[0022] The step (c) preferably includes the step of: (c1) replacing generic components generated in the previous steps with actual components selected from available devices. The generic components to be replaced include, for example, the optical amplifying units placed in the link. The step (c) preferably further includes the step of generating multiple solutions using the available devices, that can be ranked according to additional criteria. The additional criteria can be, for example, the total cost of the link. The step (c) can also include the step of determining the powers of all signal channels at the receiver, and inserting a receiver preamplifier if required to make the signals compatible with the receiver. The requirements for compatibility can include, for example, sufficient received power and sufficient received signal-to-noise ratio.

[0023] The step (d) preferably can include the step of optimising the pre-amplification gain and input channel power for the link. Further, this can include the step of: equalising the powers of input channels to obtain a substantially flat spectrum across the transmission bandwidth.

BRIEF DESCRIPTION OF THE DRAWINGS

[0024] Preferred and other embodiments of the present invention will now be described with reference to the accompanying drawings in which:

[0025]FIG. 1 illustrates schematically an example optical link that can be synthesized with the embodiment of the present invention;

[0026]FIG. 2 illustrates a flow chart of the steps in the preferred embodiment;

[0027]FIG. 3 illustrates the process ofjoining fiber segments;

[0028]FIG. 4 illustrates the process of bypassing facility locations; and

DETAILED DESCRIPTION OF THE PREFERRED AND OTHER EMBODIMENTS

[0029] In the preferred embodiment, there is disclosed a series of algorithms for determining the layout of optical components within an optical transmission link.

[0030] Initially, the embodiments will be discussed with reference to an example system of an optical link. In the example case, it is desired to form an optical link carrying multiple wavelength channels between two communication points and having a low dispersion. This kind of link is often known as an Optical Multiplex Section (OMS). A schematic illustration of the link is shown 100 in FIG. 1. A series of channels are input multiplexed 102 at one end of the link and are de-multiplexed 103 at the other end. The link includes a number of sets 105, 106 of components with each set including a portion of single mode fiber 108, a first amplifier 109, a dispersion compensating module 110 and a second amplifier 111. The dispersion compensating module may comprise a section of dispersion compensating fiber, or some other means of dispersion compensation. A further amplifier comprising a de-multiplexer preamplifier 115 is also provided.

[0031] The strategy is to place one or more sequences of these components—a Single Mode Fiber (SMF), an Amplifier (Amp), a Dispersion Compensating Module (DCM), an Amplifier (Amp)—between the Multiplexer (MUX) 102 and the De-Multiplexer (DeMUX) 103. The user is assumed to be able to specify the length and other parameters of each span of SMF e.g., 108. The Multiplexer input is connected to a Transmitter and the De-Multiplexer output is connected to a receiver.

[0032] The relevant input data can include necessary link and technology specifications including details on catalogue parts. Preferably, the system is able to generate all the signals propagated in the link and other relevant information including the location of each amplifier, its gain, the signal to noise ratio (SNR) and the transmitter power for each optical channel of the optical link.

[0033] Turning now to FIG. 2, there is illustrated a flow chart 200 of the steps of the preferred embodiment. The first step 202 is to calculate the required optical signal to noise ratio (OSNR) at the output of the link. An optimisation procedure is then used to bypass redundant facility locations by representing sections of the link using equivalent fibers 204. Next, the maximum acceptable transmission power is calculated 206. This power is checked against assumptions made previously in the process, and if these assumptions are found to be invalid the process can perform further iterations 216 using the revised value of maximum power. If the required power cannot be achieved with, e.g., the physical constraints of available devices, then the synthesis fails at this point 218. Next, the dispersion compensation components are added 28 and parameters calculated for those components. Next, a set of solutions is generating using available components to replace generic components placed in the previous steps of the algorithm 210. If no solution is found using available components, the synthesis fails at this point 220. The transmitted and received powers are then optimized 212 using transmitter pre-emphasis, and an optional receiver preamplifier, so as to produce output 214. Each of the steps 202-212 is described in more detail below.

[0034] The desired output of the algorithm is the generation of all the signals propagated in the link and other relevant information such as the location of each amplifier, the gain of each amplifier, SNR and the transmitter power for each optical channel of the optical link.

[0035] 1 Method for Calculation of Equivalent Fiber Representations (21 of FIG. 2)

[0036] An essential requirement of the Optical Multiplex Section synthesis technique of the present invention is the ability to reduce multiple concatenated actual optical fiber spans to a single effective fiber span, or equivalent fiber representation. Methods for performing this reduction are described in the following subsections.

[0037] 1.1 First Algorithm

[0038] As a first step, the following algorithm is implemented to calculate the parameters of the “effective fiber” between Facility locations. This algorithm is illustrated by way of example with reference to FIG. 3, which illustrates a series of fiber segments e.g. 31, 32 between facilities e.g., 30, 35. The algorithm can consist of the following steps:

[0039] 1. Start from the leftmost Facility location 30.

[0040] 2. Start from the left Facility location 30.

[0041] 3. Take the left-most two fiber segments e.g. 31, 32.

[0042] 4. Make a new fiber segment 33 by adding the lengths of the two fiber segments.

[0043] 5. Calculate the left-hand side variables in accordance with the definitions and equations outlined below.

[0044] 6. Use the new fiber segment 33 as the left-most fiber segment.

[0045] 7. Go to Step 3 until the length of the resulting fiber 34 reaches the next Facility 35 location.

[0046] 8. Make this Facility 35 the left-most Facility location.

[0047] 9. Go to Step 2 until there are no more Facility locations left.

[0048] 10. Verify the correctness of the results of steps 1 to 9 of the algorithm for each section with the formulas as set out in Section.

[0049] The following symbols are used to represent the actual and effective fiber sections. Symbol Description fiber_(i) i^(th) fiber fiber_(i−1) (i + 1)^(th) fiber fiber_(eff,i) Effective fiber made using fiber_(i) and fiber_(i+1)

[0050] Wherein each actual fiber has the following parameters. Symbol Description I Length A_(eff) Effective cross-section area D Dispersion coefficient DS Dispersion-slope coefficient D_(PMD) PMD coefficient α loss coefficient α_(P) loss coefficient for Raman pump band α_(PDL) PDL coefficient γ non-linear coefficient n₂ non-linear Kerr coefficient g_(Raman) Raman gain slope

[0051] The parameters of the effective fiber, fiber_(eff,i), obtained by concatenating the two fiber segments fiber_(eff,i−1) and fiber_(i) can be calculated according to the following formulae: $\begin{matrix} {{{fiber}_{{eff},\quad i}.l} = {{{fiber}_{i}.l} + {{fiber}_{{{eff},\quad i} - 1}.l}}} \\ {{{fiber}_{{eff},\quad i}.A_{eff}} = \left( \frac{{{{fiber}_{i}.A_{eff}} \cdot {{fiber}_{i}.l}} + {{{fiber}_{{{eff},\quad i} - 1}.A_{eff}} \cdot {{fiber}_{{{eff},\quad i} - 1}.l}}}{{{fiber}_{i}.l} + {{fiber}_{{{eff},\quad i} - 1}.l}} \right)} \\ {{{fiber}_{{eff`},\quad i}.D} = \left( \frac{{{{fiber}_{i}.D} \cdot {{fiber}_{i}.l}} + {{{fiber}_{{{eff},\quad i} - 1}.D} \cdot {{fiber}_{{{eff},\quad i} - 1}.l}}}{{{fiber}_{i}.l} + {{fiber}_{{{eff},\quad i} - 1}.l}} \right)} \\ {{{fiber}_{{eff},\quad i}.{DS}} = \left( \frac{{{{fiber}_{i}.{DS}} \cdot {{fiber}_{i}.l}} + {{{fiber}_{{{eff},\quad i} - 1}.{DS}} \cdot {{fiber}_{{{eff},\quad i} - 1}.l}}}{{{fiber}_{i}.l} + {{fiber}_{{{eff},\quad i} - 1}.l}} \right)} \\ {{{fiber}_{{eff},\quad i}.D_{PMD}} = \sqrt{\left( \frac{{\left( {{fiber}_{i}.D_{PMD}} \right)^{2} \cdot {{fiber}_{i}.l}} + {\left( {{fiber}_{{{eff},\quad i} - 1}.D_{PMD}} \right)^{2} \cdot {{fiber}_{{{eff},\quad i} - 1}.l}}}{{{fiber}_{i}.l} + {{fiber}_{{{eff},\quad i} - 1}.l}} \right)}} \\ {{{fiber}_{{eff},\quad i}.\alpha} = \left( \frac{{{{fiber}_{i}.\alpha} \cdot {{fiber}_{i}.l}} + {{{fiber}_{{{eff},\quad i} - 1}.\alpha} \cdot {{fiber}_{{{eff},\quad i} - 1}.l}} + {\ln ({SpliceLoss})}}{{{fiber}_{i}.l} + {{fiber}_{{{eff},\quad i} - 1}.l}} \right)} \\ {{{fiber}_{{ef},\quad {if}}.\alpha_{P}} = \left( \frac{{{{fiber}_{i}.\alpha_{P}} \cdot {{fiber}_{i}.l}} + {{{fiber}_{{{eff},\quad i} - 1}.\alpha_{P}} \cdot {{fiber}_{{{eff},\quad i} - 1}.l}} + {\ln ({SpliceLoss})}}{{{fiber}_{i}.l} + {{fiber}_{{{eff},\quad i} - 1}.l}} \right)} \\ {{{fiber}_{{eff},\quad i}.\gamma} = {\left( \frac{{fiber}_{{eff},\quad i}.\alpha}{1 - {\exp \left( {{- {{fiber}_{{eff},\quad i}.\alpha}} \cdot {{fiber}_{{eff},\quad i}.l}} \right)}} \right) \times}} \\ {{~~~}\left( {{{{fiber}_{i}.\gamma} \cdot \left( \frac{1 - {\exp \left( {{- {{fiber}_{i}.\alpha}} \cdot {{fiber}_{\quad i}.l}} \right)}}{{fiber}_{\quad i}.\alpha} \right)} + {{{fiber}_{{{eff},\quad i} - 1}.\gamma} \cdot {\exp \left( {{- {{fiber}_{i}.\alpha}} \cdot {{fiber}_{i}.l}} \right)} \cdot \left( \frac{1 - {\exp \left( {{- {{fiber}_{{{eff},\quad i} - 1}.\alpha}} \cdot {{fiber}_{{{eff},\quad i} - 1}.l}} \right)}}{{fiber}_{{{eff},\quad i} - 1}.\alpha} \right)}} \right)} \\ {{{{fiber}_{{eff},\quad i}.n_{2}} = \left( \frac{{{fiber}_{{eff},\quad i}.\gamma} \cdot c \cdot {{fiber}_{{eff},\quad i}.A_{eff}}}{2{\pi ~ \cdot f_{c}}} \right)}\quad} \\ {{{fiber}_{{eff},\quad i}.g_{Raman}} = {\frac{{{{fiber}_{{eff},\quad i}.\alpha_{P}} \cdot {fiber}_{{eff}\quad}},\quad {i.A_{eff}}}{\left( {1 - {\exp \left( {{- {{fiber}_{{eff}{,\quad}\quad i}.\alpha_{P}}} \cdot {{fiber}_{{eff},\quad i}.l}} \right)}} \right)} \times \begin{pmatrix} {{\left( \frac{1 - {\exp \left( {{- {{fiber}_{i}.\alpha_{P}}} \cdot {{fiber}_{i}.l}} \right)}}{{{fiber}_{i}.\alpha_{P}} \cdot {{fiber}_{i}.A_{eff}}} \right) \cdot {\exp \left( {{- {{fiber}_{{{eff},\quad i} - 1}.\alpha_{P}}} \cdot {{fiber}_{{{eff},\quad i} - 1}.l}} \right)} \cdot {{fiber}_{i}.g_{Raman}}} +} \\ {\left( \frac{1 - {\exp \left( {{- {{fiber}_{{{eff},\quad i} - 1}.\alpha_{P}}} \cdot {{fiber}_{{{eff},\quad i} - 1}.l}} \right)}}{{{fiber}_{{{eff},\quad i} - 1}.\alpha_{P}} \cdot {{fiber}_{{{eff},\quad i} - 1}.A_{eff}}} \right){{fiber}_{{{eff},\quad i} - 1}.g_{Raman}}} \end{pmatrix}}} \\ {{{fiber}_{{eff},\quad i}.P_{{out} - \max}} = {\min \left( {{{fiber}_{i}.P_{{out} - {\max \quad}}},\quad {{fiber}_{{{eff},\quad i} - 1}.P_{{out} - \max}}} \right)}} \end{matrix}$

[0052] 1.2 Alternative Algorithm

[0053] The first algorithm described in the previous section generates the parameters of the effective fiber by successively concatenating the effective fiber representation with additional actual fiber segments, The advantage of this approach in a specific embodiment is that a single function that concatenates two segments can be implemented, and then used whenever concatenation of two or more actual fiber segments is required. This avoids duplication of functionality within the implementation, which is a possible source of error or inconsistency. However, the successive concatenation method is mathematically equivalent to an alternative approach, in which the effective fiber representation is calculated from all segments via a set of monolithic computations. This alternative approach is described in this section.

[0054] The following symbols are now used to represent the actual fiber sections, and the overall effective fiber Symbol Description fiber_(i) i^(th) fiber from left N Total number of fiber segments fiber_(eff) Effective fiber made using N fibers

[0055] Wherein each actual fiber again has the following parameters. Symbol Description I Length A_(eff) Effective cross-section area D Dispersion coefficient DS Dispersion-slope coefficient D_(PMD) PMD coefficient α loss coefficient α_(P) loss coefficient for Raman pump band α_(PDL) PDL coefficient γ non-linear coefficient n₂ non-linear Kerr coefficient g_(Raman) Raman gain slope

[0056] Suppose the above fibers are connected serially to form the following effective fiber: Symbol Type Name Description fiber_(eff) Fiber Effective fiber made using fiber₁. . .fiber_(N)

[0057] Then the parameters of the effective fiber can be calculated directly using the following formulae: $\begin{matrix} {{{fiber}_{eff}.l} = {\sum\limits_{i = 1}^{N}\quad {{fiber}_{i}.l}}} \\ {{{fiber}_{eff}.A_{eff}} = \left( \frac{\sum\limits_{i = 1}^{N}\quad {{{fiber}_{i}.A_{eff}} \cdot {{fiber}_{i}.l}}}{\sum\limits_{i = 1}^{N}\quad {{fiber}_{i}.l}} \right)} \\ {{{fiber}_{eff}.D} = \left( \frac{\sum\limits_{i = 1}^{N}\quad {{{fiber}_{i}.D} \cdot {{fiber}_{i}.l}}}{\sum\limits_{i = 1}^{N}\quad {{fiber}_{i}.l}} \right)} \\ {{{fiber}_{eff}.{DS}} = \left( \frac{\sum\limits_{i = 1}^{N}\quad {{{fiber}_{i}.{DS}} \cdot {{fiber}_{i}.l}}}{\sum\limits_{i = 1}^{N}\quad {{fiber}_{i}.l}} \right)} \\ {{{fiber}_{eff}.D_{PMD}} = \sqrt{\left( \frac{\sum\limits_{i = 1}^{N}\quad {\left( {{fiber}_{i}.D_{PMD}} \right)^{2} \cdot {{fiber}_{i}.l}}}{\sum\limits_{i = 1}^{N}\quad {{fiber}_{i}.l}} \right)}} \\ {{{fiber}_{eff}.\alpha} = \left( \frac{\left( {\sum\limits_{i = 1}^{N}\quad {{{fiber}_{i}.\alpha} \cdot {{fiber}_{i}.l}}} \right) + {\left( {N - 1} \right) \cdot {\ln ({SpliceLoss})}}}{\sum\limits_{i = 1}^{N}\quad {{fiber}_{i}.l}} \right)} \\ {{{fiber}_{eff}.\alpha_{P}} = \left( \frac{{\sum\limits_{i = 1}^{N}\quad {{{fiber}_{i}.\alpha_{P}} \cdot {{fiber}_{i}.l}}} + {\ln ({SpliceLoss})}}{\sum\limits_{i = 1}^{N}\quad {{fiber}_{i}.l}} \right)} \\ {{{fiber}_{eff}.\gamma} = {\left( \frac{{fiber}_{eff}.\alpha}{1 - {\exp \left( {{- {{fiber}_{eff}.\alpha}} \cdot {{fiber}_{eff}.l}} \right)}} \right) \times \left( {\sum\limits_{i = 1}^{N}\quad \left( {{{fiber}_{i}.\gamma} \cdot {\exp \left( {- {\sum\limits_{j = 1}^{i - 1}\quad {{{fiber}_{j}.\alpha} \cdot {{fiber}_{j}.l}}}} \right)} \cdot \left( \frac{{1 - {\exp \left( {{- {{fiber}_{i}.\alpha}} \cdot {{fiber}_{i}.l}} \right)}}~}{{fiber}_{i}.\alpha} \right)} \right)} \right)}} \\ {{{fiber}_{eff}.n_{2}} = \left( \frac{{{fiber}_{eff}.\gamma} \cdot c \cdot {{fiber}_{eff}.A_{eff}}}{2{\pi \cdot f_{c}}} \right)} \\ {{{fiber}_{eff}.g_{Raman}} = {\frac{{{fiber}_{eff}.\alpha_{P}} \cdot {{fiber}_{eff}.A_{eff}}}{\left( {1 - {\exp \left( {{- {{fiber}_{eff}.\alpha_{P}}} \cdot {{fiber}_{eff}.l}} \right)}} \right)} \times \left( {\sum\limits_{i = 1}^{N}\quad \left( {\left( \frac{1 - {\exp \left( {{- {{fiber}_{i}.\alpha_{P}}} \cdot {{fiber}_{i}.l}} \right)}}{{{fiber}_{i}.\alpha_{P}} \cdot {{fiber}_{i}.A_{eff}}} \right) \cdot {\exp \left( {- {\sum\limits_{j = {i + 1}}^{N}\quad {{{fiber}_{j}.\alpha_{P}} \cdot {{fiber}_{j}.l}}}} \right)} \cdot {{fiber}_{i}.g_{Raman}}} \right)} \right)}} \\ {{{{fiber}_{{eff}\quad}.P_{{out} - \max}} = {\min \left( {{{fiber}_{1}.P_{{out} - \max}}{\quad,\quad}\quad \ldots \quad,\quad {{fiber}_{N}.P_{{out} - \max}}} \right)}}\quad} \end{matrix}$

[0058] 2 Link Synthesis Method

[0059] The preferred objective is to synthesize a WDM Optical Multiplex Section that satisfies a predetermined required performance objective. In this example, the objective is that, for all channels transmitted in the link, the probability of a bit being received in error is less than some user-specified bit-error rate, BER. A number of steps are used to achieve this objective, as described in the following subsections.

[0060] 2.1 Calculate required OSNR (202 of FIG. 2)

[0061] (a) Calculate Q for Given BER

[0062] The first step in the synthesis process is to convert the user-specified BER into an equivalent performance metric known as the Q-factor, or simply Q, which is better suited for use in the subsequent calculations. The Q is computed from the BER using the following equation: ${BER} = {\frac{1}{\sqrt{2\pi}Q}{\exp \left( {- \frac{Q^{2}}{2}} \right)}}$

[0063] Note that Q is higher for a channel with better performance (i.e. lower BER).

[0064] (b) Adjust Q to Account for Four-Wave Mixing (FWM) Penalty

[0065] When multiple WDM channels are transmitted at high power levels through a single optical fiber, a nonlinear process known as four-wave mixing (FWM) results in crosstalk between the channels. This causes a penalty at the receiver, whereby the BER will be higher than that predicted if the effect of FWM is not taken into account. In order to ensure that this penalty does not result in a failure to meet the BER objective, the value of Q found in the previous step is preferably adjusted upwards to provide a margin for the FWM penalty.

[0066] The adjusted value of Q is determined using the following set of procedures.

[0067] An initial procedure is used to define the system properties:

[0068] 1. Consider M spans of fiber with the following properties:

[0069] γj=Nonlinear index of j^(th) fiber span in units of W⁻¹m⁻¹

[0070] Dj=Dispersion ofj^(th) fiber span in units of s m⁻²

[0071] 2. Consider N optical signals, with equal power, equally spaced on a frequency grid, with the following properties:

[0072] λ=Center wavelength

[0073] Δf=Frequency spacing of channels

[0074] Pj =Power at start of each fibre section

[0075] 3. For the system, compute:

σ₁=1.4×(1−1.7N^(−0.75))

[0076] 4. For each section of fiber compute: $\begin{matrix} {k_{j} = \frac{P_{j}\gamma \quad c}{\left. {\pi\lambda}^{2} \middle| D_{j} \middle| {\Delta \quad f^{2}} \right.}} \\ {{F_{j}(p)} = {- \frac{2}{{k_{j}^{2}\sigma_{1}^{2}} + 2}}} \end{matrix}$

[0077] 5. For the entire span compute:

[0078] ${F(p)} = {\prod\limits_{j}^{\quad}{F_{j}\quad (p)}}$

ƒ(S)=ℑ⁻¹ {F(p)}

[0079] where ℑ⁻¹ represents the Inverse Fourier Transform and ƒ(S) is the probability density function (PDF) of the parameter S.

[0080] 6. Define p to be the ratio of noise during marks and spaces for the undistorted system, ρ=σ_(mark,0)/σ_(space,0). ρ typically takes on a value of approximately 2.

[0081] 7. ${{Define}\quad {{erfc}(Q)}} = {\frac{1}{\sqrt{\pi}}{\int_{Q/\sqrt{2}}^{\infty}{^{- \frac{^{2}}{2}}\quad {t}}}}$

[0082] Then the following procedure defines a function Q_(avg)=FWMQpenalty(Q₀):

[0083] 1. Given a Q value for the undistorted system (without FWM considered), Q₀, a BER value can be calculated as follows taking FWM into account: ${{BER}\left( {S,\quad \rho,\quad Q_{0}} \right)} = {\frac{1}{2}\left\lbrack {{{erfc}\left( {Q_{0} \times \frac{{\left( {1 - S} \right)\left( {\rho + 1} \right)} - 1}{\sqrt{1 + {\left( {\rho^{2} - 1} \right)\left( {1 - S} \right)}}}} \right)} + {{erfc}\left( Q_{0} \right)}} \right\rbrack}$

[0084] 2. The average BER for the system with FWM considered is computed, weighted by the PDF of S, ƒ(S).

BER _(avg) =∫BER(S,ρ,Q ₀)ƒ(S)dS

[0085] 3. An average Q value is computed by solving the following equation for Q_(avg):

BER _(avg) =erƒc(Q_(avg)).

[0086] A solution to the equation Q_(desired)=FWMQpenalty(Q₀) is required for Q₀ε[4,12] using a specified desired value of Q=Qdesired, being equal to the initial Q-factor determined from the user-specified BER. The adjusted value Q₀, if it exists, will be greater than the target Q_(desired), thus providing a margin that will allow for the degradation in Q-factor due to FWM. The required solution is obtained using the following procedure:

[0087] 1. The function FMWQpenalty is evaluated at values Q₀=4 and Q₀=12,

Q _(test,1) =FMWQpenalty(4)

Q _(test,2) =FWMQpenalty(12)

[0088] 2. If the desired value does not satisfy Q_(test,1)<Q_(desired<Q) _(test,2), then no solution can be obtained for the specified fibre and signal parameters. A lower value of k is needed to obtain the desired Q value. Since k is proportional to the signal power level, a solution may be obtained if it is acceptable to the user to use a lower power level P for all channels.

[0089] 3. Otherwise, if the desired value does satisfy Q_(test,1)<Q_(desired)<Q_(test,2), then a solution for Q₀ is found in this range using the Newton-Raphson numerical rootfinding method.

[0090] The resulting value of Q₀ is the FWM-adjusted Q value that is used in the following steps of the link synthesis procedure.

[0091] (c) Calculate Required OSNR at Receiver

[0092] The optical signal-to-noise ratio (OSNR) is the ratio of the optical signal power to the accumulated optical noise (added in optical amplifiers) within the optical bandwidth of a signal channel, as measured at the receiver. The required OSNR at the receiver is determined from the adjusted Q value obtained at the previous step by solving the following equation: $Q \approx \frac{2 \times {OSNR} \times \left( \sqrt{\alpha} \right)}{\sqrt{\left( {1 - \frac{1}{2\alpha}} \right)} + \sqrt{\left( {1 - \frac{1}{2\alpha}} \right) + {4 \times {OSNR}}}}$

[0093] where α is a user-specified parameter representing the ratio of the optical bandwidth of each signal channel to the bit-rate of the data carried on the channel.

[0094] There are a number of additional processes and factors that can further degrade the quality of the signal at the receiver. Each of these effects can be accounted for empirically by assignment of an effective OSNR margin. The complete set of margins constitutes a heuristic method to include all sources of degradation. The margins are preferably specified via a corresponding set of user-supplied variables such as those listed in the following table. Variable Name Description heuTxChrpMargin Margin to account for excess bandwidth of the transmitter compared with a chirpless modulator heuTxElecMargin Margin to account for imperfect (i.e. not flat amplitude, not linear phase) response of the transmitter drive circuitry and laser parasitics heuRxElecMargin Margin to account for imperfect amplitude and phase response of the electrical receiver heuRxSenTiltMargin Margin to account for signal reduction in wavelength-dependence of receiver heuTxAgeMargin Margin to account for transmitter aging heuRxAgeMargin Margin to account for receiver aging heuMuxBwPhMargin Margin to account for imperfect (not flat) amplitude and phase (non-linear) response of optical multiplexer on desired channel heuMuxPolMargin Margin to account for polarization dependence of multiplexer heuMuxAgeMargin Margin to account for degradation of multiplexer response due to aging heuMuxXtalkMargin Margin to account for crosstalk from adjacent channels in an imperfect multiplexer heuDeMuxBwPhMargin Margin to account for imperfect (not flat) amplitude and phase (non-linear) response of optical de-multiplexer on desired channel heuDeMuxPolMargin Margin to account for polarization dependence of de-multiplexer heuDeMuxAgeMargin Margin to account for degradation of de-mul- tiplexer response due to aging heuDeMuxXtalkMargin Margin to account for crosstalk from adjacent channels in an imperfect de-multiplexer heuFiberRepairMargin Margin to account for loss due to repairs of fiber in the optical link heuFiberAgeMargin Margin to account for increase of fiber loss due to aging heuMonitorMargin Margin to account for power loss due to per- formance monitoring heuDispMargin Margin to account for linear dispersion in the whole fiber plant heuDispCompMargin Margin to account for signal degradation caused by dispersion-compensating equiptment with imperfect characteristics heuSPMMargin Margin to account for SPM induced effects heuRamanMargin Margin to account for power transfer between channels by Raman Scattering in the WDM system heuBrillouinMargin Margin to account for power loss by Brillouin scattering in the WDM system heuPMDMargin Margin to account for PMD heuPMLMargin Margin to account for PML heuAmpGTiltMargin Margin to account for signal reduction in worst channels due to gain tilt heuAmpPolGMargin Margin to account for polarization depend- ent gain of amplifiers heuAmpPolHoleMargin Margin to account for polarization hole burning in amplifiers heuAmpAgeMargin Margin to account for gain reduction in am- plifiers due to aging heuFECBERGain BER improvement due to Forward Error Correction (FEC)

[0095] An adjusted OSNR, OSNR_(adj), can then be calculated from the unadjusted OSNR using the following formulae:

[0096] FECSNRGain=2_(×)ln(heuFECBERGain) $\begin{matrix} {{SNRMargin} = {{heuTxChirpMargin} \times}} \\ {{{heuTxElecMargin} \times}} \\ {{{heuRxElecMargin} \times}} \\ {{{heuRxSenTiltMargin} \times}} \\ {{{heuTxAgeMargin} \times}} \\ {{{heuRxAgeMargin} \times}} \\ {{{heuMuxBwPhMargin} \times}} \\ {{{heuMuxPolMargin} \times}} \\ {{{heuMuxAgeMargin} \times}} \\ {{{heuMuxXtalkMargin} \times}} \\ {{{heuDeMuxBwPhMargin} \times {heuDeMuxPolMargin} \times}} \\ {{{heuDeMuxAgeMargin} \times}} \\ {{{heuDeMuxXtalkMargin} \times}} \\ {{{heuFiberRepairMargin} \times}} \\ {{{heuFiberAgeMargin} \times}} \\ {{{heuMonitorMargin} \times}} \\ {{{heuDispMargin} \times}} \\ {{{heuDispCompMargin} \times}} \\ {{{heuSPMMargin} \times}} \\ {{{heuRamanMargin} \times}} \\ {{{heuBrillouinMargin} \times}} \\ {{{heuPMDMargin} \times}} \\ {{{heuPMLMargin} \times}} \\ {{{heuAmpGTiltMargin} \times}} \\ {{{heuAmpPolGMargin} \times}} \\ {{{heuAmpPolHoleMargin} \times}} \\ {{{heuAmpAgeMargin}/}} \\ {{FECSNRGain}} \end{matrix}$

[0097] OSNR_(adj)=OSNR_(×)SNRMargin

[0098] 2.2 Bypass Redundant Facility Locations Using Equivalent Fiber Representation (204 of FIG. 2)

[0099] Next an optimization algorithm is used to bypass redundant facility locations. The object of the optimization is to eliminate the largest number of facility locations that can be found, that are unnecessary in order to meet the link design objectives. Advantageously, each facility location that can be eliminated in this way represents a saving in equipment and maintenance expenses that substantially reduces the installation and operating costs of the link.

[0100] (a) Constraints

[0101] A facility location can be eliminated if, after concatenating the fiber sections adjacent to the facility, a number of constraints are satisfied. These constraints can be expressed as follows.

[0102] First, a number of fixed design constraints and objectives must be considered.

[0103] there is a maximum gain G_(max) achievable by any optical amplifier that can be placed at a facility location, and thus the total loss between any two facility locations must not exceed this value;

[0104] the user may wish to limit the maximum dispersion to some fraction of the theoretical dispersion length, which can be achieved by setting a heuristic “dispersion length contraction factor” represented by a variable heuKi;

[0105] the user may wish to limit the absolute maximum span length, which can be achieved by setting a heuristic control represented by a variable heuMaxSpanLength;

[0106] there is a maximum acceptable level of total optical noise that can be contributed by all amplifiers combined as determined by the required OSNR, which can be achieved by constraining the total gain in the link to a fixed maximum that depends upon the noise figure NF of the optical amplifiers, and that is represented by a variable MaxTotalGain;

[0107] there is a target total link capacity (in bits per second), represented by a variable linkCapacity, which it is the objective of the design process to achieve;

[0108] there is a technology limit on the spacing between adjacent wavelength channels, represented by a variable techChannelSpacing.

[0109] Next, it is a requirement that the variation in loss (or, equivalently, in the compensating gain required) for each segment is minimum, as well as minimizing the total number of segments (which is assumed to be an overall objective). In other words, it is required to find a particular configuration in which the lengths of all segments are as uniformly distributed as possible. Advantageously, a design that meets these requirements will be easier to deploy and to maintain.

[0110] The above requirements can be expressed mathematically, in a form suitable for implementation in the preferred embodiment as described in the following section. The facility bypass algorithm constructs sets of concatenated fiber segments that can be expressed using the following equivalent set equation:

[0111] ψ_(i)=[n_(1i), n_(1i)+1,n_(1i)+2, . . . ,n_(mi)]with 1≦m≦N ,1≦i≦N, and n_(1j)ε[1,N] are consecutive integers for a given ψ_(i). These sets have the property that: ${{\psi_{i}\bigcap\limits_{\forall{i \neq j}}\psi_{j}} = {{\varnothing \quad {and}\quad \underset{i}{U}\psi_{i}} = \left\lbrack {1,2,\ldots \quad,N} \right\rbrack}},$

[0112] i.e. all the fiber segments from the previous phase are used, and the set is ordered in the sense that if i ≦j then ψ_(i)≦ψ_(j):ψ_(i)≠,ψ_(j)≠ where relational operators are operated pair-wise between the sets. The fixed design constraints and objectives described previously can then be expressed in the form:

G_(i)≦G_(max)

[0113] where $G_{i} = {\exp \left( {\sum\limits_{j \in \psi_{i}}^{\quad}\quad {{\overset{\_}{a}}_{j}{\overset{\_}{l}}_{j}}} \right)}$

[0114] is the gain required to compensate for the loss in the i^(th) concatenated segment, that is provided by the line amplifiers at the end of each segment, e.g. 109 of FIG. 1. $\begin{matrix} {{\sum\limits_{j \in \psi_{i}}^{\quad}\quad {\overset{\_}{l}}_{j}} \leq {{heuKi} \times {\overset{\_}{l}}_{{UB}\quad \psi_{\quad_{i}}}}} \\ {{\sum\limits_{j \in \psi_{i}}^{\quad}\quad {\overset{\_}{l}}_{j}} \leq {heuMaxSpanLength}} \\ {{\prod\limits_{i}^{\quad}\quad G_{i}} \leq {MaxTotalGain}} \end{matrix}$

[0115] while the minimization in variation of lengths can be achieved by minimizing the following quantity: V = ( ∑ i = 1     ( L ψ i - ∑ i = 1     L ψ i / ) 2 ) /

[0116] In the above expressions, the following definitions of symbols hold:

[0117] L_(ψ) _(i) is the total loss of the fiber segment in section ψ_(i;)

[0118]

is the number of sections;

[0119] {overscore (l)}_(UBψdi i)={square root}{square root over (l)}_(Dψ) _(i) ×{overscore (l)}_(NLψ) _(i) is the length upper bound; $\begin{matrix} {{{\overset{\_}{l}}_{{NL}_{\psi_{i}}} = {\frac{1}{{\overset{\_}{\gamma}}_{\psi_{i}}{\overset{\_}{P}}_{\psi_{i}}^{thresh}}\left( \frac{linkCapacity}{2*{techChannelSpacing}} \right)\quad {is}\quad {the}\quad {``{{nonlinear}\quad {length}}"}}}\quad;} \\ {{{\overset{\_}{l}}_{D_{\psi_{i}}} = {\frac{2\pi \quad c \times {heuTo}^{2}}{\left. \lambda_{ref}^{2} \middle| {\overset{\_}{D}}_{\psi_{i}} \right|}\quad {is}\quad {the}\quad {dispersion}\quad {length}}};} \end{matrix}$

[0120] The parameter {overscore (P)}₁₀₄ _(i) ^(thresh) that appears in the expression for the nonlinear length is the maximum optical power that must be handled by the link, which in the preferred embodiment is set to the value of P_(out-max) determined in the equivalent fiber representation calculation. The parameter heuTo that appears in the expression for the dispersion length is a heuristic parameter that represents an equivalent pulse width (full width half maximum) for the transmitted data bits. This parameter is determined according to the characteristics of the transmitter(s).

[0121] (b) Divide Link into Independent Blocks

[0122]FIG. 4 illustrates a link comprising a series of fiber segments e.g. 41, 42 located between facilities e.g. 40, 45. The first step in the redundant facility elimination algorithm is to divide the link into independent sub-links. A pair of independent sub-links can be identified as follows:

[0123] 1. Starting from the left-most two adjacent fiber sections e.g. 41, 42.

[0124] 2. Determine the total loss of the effective fiber obtained by concatenating the two sections.

[0125] 3. If the total loss thus obtained exceeds the maximum gain G_(max) that can be achieved at a single facility location, then it will clearly be impossible for any single concatenated section to contain the two original fiber sections. In this case, the facility location between the two fiber section will be required, and it is possible to consider the two sub-links comprising all sections to the left of the facility location, and all sections to the right of the facility location, independently.

[0126] 4. If the right-hand end of the link has not been reached, return to step 2, considering the concatenation of the second fiber section from this step, e.g. 42, with the next adjacent section to its right.

[0127] (c) Eliminate Redundant Facility Locations

[0128] Redundant facility locations are eliminated on each sub-link according to the following algorithm, again with reference to FIG. 4:

[0129] 1. Start from the left-most facility location e.g. 40.

[0130] 2. Take the left-most two fiber sections e.g. 41, 42.

[0131] 3. Make a new fiber section e.g. 43 by concatenating the two fiber sections.

[0132] 4. If the properties of the new fiber section 43 satisfy the constraints set out in (a) above, then make this new fiber section the left-most fiber section and go to step 2.

[0133] 5. Take the next left-most facility location 45 and go to step 2.

[0134] 2.3 Calculate Maximum Transmission Power Allowing for Nonlinear Effects (206 of FIG. 2)

[0135] (a) Calculate Maximum Allowed Transmission Power per Channel at Transmitter Considering FWM and ASE

[0136] In this step, it is necessary to verify that the design synthesized by the previous steps is viable. The combination of FWM and optical noise added by the amplifiers places an upper limit on the power in the link, above which the assumptions made previously regarding the impact of nonlinearities are invalid. This upper limit P_(FWM) can be calculated as follows: $P_{FWM} = \sqrt{1.5\quad {hvB}_{0}{{NF}\left( {\sum\limits_{i = 1}^{N}\quad G_{i}} \right)}\sqrt{\frac{SNR}{{\sum\limits_{i = 1}^{N}\quad \frac{k_{i}^{2}k_{i}^{2}}{G_{i}}}}}}$

[0137] where

[0138] h=6.626176 ×10⁻³⁴Js

[0139] c=299792458 m/s

[0140] NF is the noise figure set for the generic variable power amplifiers

[0141] SNR is the required signal to noise ratio of the link

[0142] {l_(i):iε[1,N]} is the length of fibers in the fiber-plan

[0143] {G_(i):iε[1,N]} is the gain of the amplifiers for compensating losses $k_{i} = \frac{\gamma_{i}c}{2{\pi \left( {{\overset{\_}{D}}_{i}l_{i}} \right)}\lambda^{2}\Delta \quad f^{2}}$

[0144] {overscore (D)}_(i) is the mean dispersion seen by the channel in the i^(th) link

[0145] λ is the mean operating wavelength: normally 1550 nm

[0146] Δƒ is the minimum channel spacing of the link

[0147] γI is the nonlinear coefficient of the i^(th) fiber $= {{\max\limits_{n \in {\lbrack{1,\quad M}\rbrack}}{\sum\limits_{\quad \underset{\underset{{r,\quad p{,\quad}\quad q} \in {\lbrack{1,\quad M}\rbrack}}{{n,\quad r} \neq {p,\quad q}}}{n = {p + q - r}}}^{\quad}\quad \frac{1}{16\left( {n - p} \right)^{2}\left( {n - q} \right)^{2}}}} \approx {{0.4 \times \left( {1 - \frac{1}{M}} \right)} - {\frac{0.46}{\sqrt{M}} \times {\exp \left( {- \frac{M}{10}} \right)}}}}$

[0148] M is the number of channels

[0149] B₀=0.7 ×α× avgBitRate

[0150] The parameter avgBitRate is the average of the individual transmitter bit-rates set out in the link specification provided by the user.

[0151] If the calculated maximum allowed power P_(FWM) exceeds the power launched into any fiber section (the maximum of which is given by P_(out-max)), then it is possible to go back and repeat the synthesis process using the value of P_(FWM) as the maximum allowed launch power.

[0152] (b) Calculate Maximum Allowed Transmission Power per Channel at Transmitter Considering SRS and ASE

[0153] The nonlinear process of Stimulated Raman Scattering (SRS) can also result in an unacceptable degradation of signal quality at the receiver, if the maximum power launched into any fiber segment exceeds a threshold power P_(SRS). SRS was not included in the synthesis process, on the assumption that this threshold power level is not exceeded. It is desirable to verify this assumption. The SRS threshold power can be calculated using the following equation: $\left. {P_{SRS} = \sqrt{\left( \frac{1.5{SNR} \times \chi \times {hvB}_{0}{{NF}\left( {\sum\limits_{i = 1}^{N}\quad G_{i}} \right)}}{\frac{1}{4}{M\left( {M - 1} \right)}\Delta \quad {f\left( {\sum\limits_{i = 1}^{N}\quad {\left( \frac{g_{Roman}}{A_{eff}} \right)_{i}\frac{l_{i}}{\ln \left( G_{i} \right)}\left( {1 - \frac{1}{G_{i}}} \right)}} \right)}} \right.}} \right)$

[0154] where:

[0155] χ is a user-specified parameter defining the acceptable SRS crosstalk level (typically 1 dB), and all other parameters have been defined previously.

[0156] If the calculated maximum allowed power P_(SRS) exceeds the power launched into any fiber section (the maximum of which is given by P_(out-max)), then it is possible to go back and repeat the synthesis process using the value of P_(SRS) as the maximum allowed launch power.

[0157] (c) Verify that Per-Channel-Output-Power is Valid The output power at each transmitter is set according to:

P_(set)=min(P_(FWM),P_(SRS),P_(Tx−max))

[0158] where P_(Tx-max) is the maximum output power that all transmitters are capable of generating.

[0159] This output power setting is then validated against the target SNR, OSNR_(adj), by verifying that the following relation holds: $\frac{P_{set}}{1.5{hvB}_{0}{{NF}\left( {\sum\limits_{\quad i}^{\quad}\quad G_{i}} \right)}} \geq {OSNR}_{adj}$

[0160] The factor of 1.5 in the denominator accounts for the additional noise introduced when dispersion compensation, and associated amplifiers, are added to the design in the next step.

[0161] If this relation does not hold, then the available output power from the user-specified transmitters is not sufficient to achieve the design objectives, and the synthesis has failed.

[0162] 2.4 Apply Dispersion Compensation (208 of FIG. 2)

[0163] (a) Determine Residual Dispersion

[0164] When high power signals are launched into an optical fiber, the nonlinear process of self-phase modulation (SPM) causes spectral distortion of the signals. Under these conditions, signal quality at the receiver can generally be improved by under-compensating for dispersion, and thus allowing a residual dispersion to exist over the link. The optimum amount of residual dispersion can be calculated using the following procedure.

[0165] 1. The SPM limited per-channel power-distance product is determined for a given set of equipment. In general it can be difficult to obtain this value by analytical means, and therefore a physical layer simulation software tool such as VPItransmissionMaker™ WDM available from VPlsystems Inc may be used. This value can be denoted by the name PowXDistProd.

[0166] 2. The residual dispersion value at the receiver can then be optimised by finding the residual dispersion for which the Q value at the receiver is minimum. Again, if no more efficient means is available, a physical layer simulation software tool such as VPItransmissionMaker™ WDM may be used. This value can be denoted by the name ResDispVal.

[0167] 3. The residual dispersion δDL_(i) at the end of each individual fiber span i of length l_(i) given the input power P per channel is then determined using the following formula: ${\delta \quad {DL}_{i}} = {- \left( \frac{{ResDispVal} \times l_{i} \times P}{{UPerImFactor} \times {PowXDistProd}} \right)}$

[0168]  where UPerImFactor is a user-supplied estimated performance improvement factor due to under-compensation, which typically assumes a value around 2.

[0169] (b) Apply Dispersion Compensation by Selection from Available Dispersion Compensating Modules

[0170] Compensation for dispersion in each fiber span is provided by dispersion compensating modules (DCM) installed at the end of the span, e.g. 110 of FIG. 1. The amount of dispersion compensation (DL)comp,J required after span i to achieve the target residual dispersion δDL_(i) is calculated according to:

(DL _(comp,i))=−D _(i) l _(i) +δDL _(i)

[0171] The objective is to select a dispersion compensating solution from a user-supplied library of available DCM components, each of which is characterized by at least its total dispersion value DC, and its loss L_(DC). The first DCM in the library satisfying the following two criteria is selected:

[0172] sgn(DC)=−sgn(D_(i)) i.e. the DCM must have opposite dispersion to the fiber span; and |DC−D_(i)l_(i)|≦ε

[0173] where ε is a user-specified fitting tolerance that defines how accurately the actual residual dispersion must approximate the target residual dispersion.

[0174] If there is no single DCM that satisfies the two criteria above, then an attempt may be made to attain the required compensation by combining DCMs with smaller total dispersion than required. If this fails, then a generic DCM may be assumed with the desired value, and with a loss L_(DC) equal to 75% of the total loss of the preceding fiber span. The properties of the generic DCM approximate the typical properties of dispersion compensating fiber (DCF), and thus the generic DCM represents a conventional dispersion compensating module constructed from a suitable length of DCF.

[0175] (c) Compensate for Loss of DCMs The loss of each DCM is compensated by inserting amplifiers with gain G_(DCM,i) equal to the loss L_(DC,i) after the DCM at the end of each fiber span i, e.g. 111 of FIG. 1.

[0176] 2.5 Generate Solutions (210 of FIG. 2)

[0177] (a) Fit Amplifiers to Generic Line Amplifiers

[0178] At this stage of the synthesis process, the link design is complete except that the line amplifiers e.g. 109, 111 of FIG. 1, which compensate for loss in each fiber span, e.g. 108 of FIG. 1, and for each DCM, e.g. 110 of FIG. 1, are generic units. The required gain at each amplifier is known, however a specific available component must be substituted for the generic device.

[0179] Available amplifiers are preferably selected from a user-supplied library of amplifier and attenuator components. The required gain at each generic amplifier is achieved by selecting one or more amplifier components followed by zero or more attenuator components from this library. Although the total gain is the same if the attenuator(s) are placed prior to the amplifiers, this configuration is not used because the degradation of signal-to-noise ratio will be higher in this case.

[0180] The maximum number of amplifiers or attenuators that can be cascaded in order to achieve the generic gain target may be limited to a user-defined maximum.

[0181] If no combination of components can be found in the library that satisfy the generic gain requirements, then the design objectives cannot be achieved, and the synthesis has failed.

[0182] (b) Generate Solutions According to User Criteria

[0183] When fitting library components (e.g. DCMs, amplifiers and attenuators) to generic requirements, there may be more than one component, or combination of components, that meet the requirements. In this case, multiple solutions may be generated. The user may be able to select from these solutions either manually, or automatically by specifying additional criteria for the design, assuming that suitable information is available in the component library. For example, a common requirement might be that a lowest cost solution is required. If cost information for each component is available in the library, then the lowest cost solution can easily be determined from multiple available solutions by calculating the total cost of each and retrieving the solution for which this total is a minimum.

[0184] 2.6 Optimize Transmitted and Received Optical Power (212 of FIG. 2)

[0185] (a) Insert Receiver Preamplifier (if required)

[0186] For each solution, the actual received average power of every channel can now be calculated using the characteristics of the components chosen from the component libraries. Due to the wavelength-dependence of amplifier gain, and component losses, the received average power may be different for different channels.

[0187] If there is insufficient received average power at any receiver, a receiver preamplifier, 115 of FIG. 1, can be introduced to increase the received power.

[0188] If the received power at any receiver is smaller than the receiver sensitivity level, then the average power is raised using a pre-amplifier until this condition is satisfied for the worst affected channel.

[0189] If the minimum SNR value at any receiver is smaller than the required SNR value, the average power is raised using a pre-amplifier until this condition is satisfied for the minimum SNR channel.

[0190] (b) Equalize Output Power

[0191] As mentioned above, once real components are introduced it is expected that the received power will be different for each wavelength channel. Variations in powers P_(i) ^(out) at the output of the link can be compensated by adjusting the corresponding transmitted powers P_(i) ^(in) at the input to the link. This technique is known as pre-emphasis.

[0192] The following algorithm is applied to compute a suitable set of pre-emphasized inputs for M channels, assuming that the desired maximum variation in received power is represented by a user-specified variable EqPowTol.

[0193] 1. Start with a set of equal initial input powers (P_(i) ^(in))_(old), and compute the corresponding output powers (P_(i) ^(out))_(old)

[0194] 2. Compute a measure of the mean-squared variation in received power for the “old” set of input and output powers according to: $\sqrt{\frac{{\sum\limits_{i = 1}^{N}\left( {\left( P_{i}^{out} \right)_{old} - {\sum\limits_{j = 1}^{N}\quad {\left( P_{j}^{out} \right)_{old}/M}}} \right)^{2}}\quad}{M}}$

[0195] 3. If the resulting measure of variation is less than EqPowTol, then a solution has been found and the algorithm terminates.

[0196] 4. Otherwise, calculate a “new” set of input powers (p_(i) ^(in))_(new) according to: $\left( P_{i}^{in} \right)_{new} = {\left( {\sum\limits_{j = 1}^{M}\quad \left( P_{j}^{in} \right)_{old}} \right)\left( \frac{\left( P_{i}^{in} \right)_{old}/\left( P_{i}^{out} \right)_{old}}{\sum\limits_{j = 1}^{M}\quad {\left( P_{j}^{in} \right)_{old}/\left( P_{j}^{out} \right)_{old}}} \right)}$

[0197] 5. If any of the resulting input powers are outside the acceptable range defined by the maximum and minimum transmitter power, and the nonlinear limits P_(FWM) and P_(SRS), then the power shall be set to the corresponding limiting value.

[0198] 6. Compute the corresponding set of output powers (P_(i) ^(out))_(new).

[0199] 7. The “new” values, now become the “old” values for the next iteration, commencing again at step 2.

[0200] 8. If this process fails to achieve a solution after a user-specified maximum number of iterations, then it shall be assumed that the target EqPowTol is not achievable.

[0201] (This may or may not imply a failure of the synthesis process, depending upon the criticality of attaining EqPowTol.)

[0202] The aforementioned example describes a method for implementing a system for minimizing the complexity of a predetermined optical arrangement. Obviously, depending on requirements, the method of the preferred embodiment can be readily extended to other arrangements.

[0203] It will be understood that the invention disclosed and defined herein extends to all alternative combinations of two or more of the individual features mentioned or evident from the text or drawings. All of these different combinations constitute various alternative aspects of the invention.

[0204] The foregoing describes embodiments of the present invention and modifications, obvious to those skilled in the art can be made thereto, without departing from the scope of the present invention. 

We claim:
 1. A method of calculating an equivalent optical fiber representation of a multiplicity of optical fibers, the method including the step of: (a) successively combining the parameters of a pair of adjacent optical fibers to compute the parameters of a combined optical fiber that has substantially the same transmission properties as the two adjacent optical fibers; wherein one of the pair of adjacent optical fibers is a previously computed combined optical fiber.
 2. A method of calculating an equivalent optical fiber representation of a multiplicity of optical fibers, the method including the step of: (a) combining the parameters of a multiplicity of optical fibers comprising an arbitrary combination of actual fibers and equivalent fiber representations, to compute the parameters of a single replacement fiber that has substantially the same transmission properties as the multiplicity of fibers.
 3. A method for optimizing the topology of an optical data transmission link, the link comprising a series of optical components interconnected by optical fiber sections, the method comprising the steps of: (a) calculating an equivalent fiber representation of the link; and (b) utilising said equivalent fiber representation to determine if any of the series of optical components can be eliminated from the link.
 4. A method as claimed in claim 3 wherein said step (b) further comprises the step of: (b1) calculating a measure of signal quality required at a receiver to satisfy a link performance criteria.
 5. A method as claimed in claim 3 wherein said step (b) further comprises the step of: (b2) calculating a maximum acceptable transmission power for said link.
 6. A method as claimed in claim 3 wherein said step (b) further comprises: (b3) determining suitable configuration values for predetermined components in the link.
 7. A method as claimed in claim 6 wherein said configuration values include at least one of the amplification factor of an optical amplifier or the total dispersion of a dispersion compensating module.
 8. A method as claimed in claim 3 further comprising the step of: (c) generating a series of components for use in the transmission link from a set of available components.
 9. A method as claimed in claim 3 further comprising the step of: (d) utilising an optimisation procedure to determine the input channel power of the optical data transmission link.
 10. A method as claimed in claim 1 wherein the step (a) further comprises the step of:(a1) progressively joining segments together between facilities and determining whether an equivalent fiber representation for the joined segment falls within first predetermined constraints.
 11. A method as claimed in claim 1 where the step (a) further comprises the step of: (a2) progressively joining segments together across facility locations and determining whether an equivalent fiber representation for the joined segment falls within predetermined constraints.
 12. A method as claimed in claim 4 wherein said step (b1) comprises the step of calculating a numerical measure of signal quality corresponding to a specified link performance.
 13. A method as claimed in claim 12 wherein the numerical measure of link performance is the probability of transmission errors occurring on the link.
 14. A method as claimed in claim 12 wherein said numerical measure comprises a Q-factor.
 15. A method as claimed in claim 4 wherein said step (b1) further includes adjusting the signal quality measure to account for additional sources of degradation in transmission.
 16. A method as claimed in claim 15 wherein said additional sources include nonlinear transmission processes such as four-wave mixing or stimulated Raman scattering.
 17. A method as claimed in claim 4 wherein said step (b1) comprises calculating an optical signal-to-noise ratio (OSNR).
 18. A method as claimed in claim 5 wherein said step (b2) further comprises calculating the maximum transmission power before which unacceptable signal degradation occurs due to nonlinear transmission effects and additive noise.
 19. A method as claimed in claim 18 wherein said additive noise comprises amplified spontaneous emission noise introduced in optical amplifiers.
 20. A method as claimed in claim 18 wherein said additive noise comprises four wave mixing or stimulated Raman scattering.
 21. A method as claimed in claim 5 further comprising the step of verifying that the output power is valid, within the operating ranges of components in the link.
 22. A method as claimed in claim 21 wherein said verification step includes verifying that the transmitters are capable of generating the required output power.
 23. A method as claimed in claim 6 further comprising the step of determining the amount of residual dispersion required to optimise the signal quality at the receiver, considering the impact of nonlinear propagation effects.
 24. A method as claimed in claim 23 wherein said nonlinear propagation effects include the self phase modulation of each transmitted channel.
 25. A method as claimed in claim 6 further comprising the step of applying dispersion compensation to achieve the required residual dispersion by selecting from available dispersion compensating modules.
 26. A method as claimed in claim 25 wherein said dispersion compensating modules can include lengths of dispersion compensating fiber.
 27. A method as claimed in claim 25 further comprising the step of compensating for the insertion loss of the dispersion compensating modules by including further optical amplifying units in the link.
 28. A method as claimed in claim 8 further comprising the step of generating multiple solutions using the available devices and ranking those solutions according to an additional predetermined criteria.
 29. A method as claimed in claim 28 wherein said additional predetermined criteria is the total cost of the link.
 30. A method as claimed in claim 9 further including the step of determining the powers of all signal channels at the receiver, and inserting a receiver preamplifier if required to make the signals compatible with the receiver.
 31. A method as claimed in claim 30 wherein said compatibility can include sufficient received power and sufficient received signal-to-noise ratio.
 32. A method as claimed in claim 9 further comprising the step of optimising the pre-amplification gain and input channel power for the link.
 33. A method as claimed in claim 32 further comprising the step of equalising the powers of input channels to obtain a substantially flat spectrum across the transmission bandwidth.
 34. A method of optimising the componentry of an optical data transmission link, the link comprising a series of optical components interconnected by a fiber network, the method comprising the steps of: (a) calculating an equivalent fiber representation for the link; (b) utilising an optimisation algorithm to eliminate redundant optical components from the link; (c) utilising an optimisation algorithm to determine the input channel power of said link. 